Einstein static Universe in hybrid metric-Palatini gravity 
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Hybrid metric-Palatini gravity is a recent and novel approach to modified theories of gravity, 
which consists of adding to the metric Einstein-Hilbert Lagrangian an f{TZ) term constructed a 
la Palatini. It was shown that the theory passes local tests even if the scalar field is very light, 
and thus implies the existence of a long-range scalar field, which is able to modify the dynamics in 
galactic and cosmological scales, but leaves the Solar System unaffected. In this work, motivated 
by the possibility that the Universe may have started out in an asymptotically Einstein static 
state in the inflationary universe context, we analyse the stability of the Einstein static Universe by 
considering linear homogeneous perturbations in the respective dynamically equivalent scalar-tensor 
representation of hybrid metric-Palatini gravity. Considering linear homogeneous perturbations, the 
stability regions of the Einstein static universe are parametrized by the first and second derivatives 
of the scalar potential, and it is explicitly shown that a large class of stable solutions exists in the 
respective parameter space, in the context of hybrid metric-Palatini gravity. 

PACS numbers: 04.20.Jb, 04.50.Kd 



I. INTRODUCTION 

The standard model of cosmology is remarkably suc- 
cessful in accounting for the observed features of the Uni- 
verse and during the last two decades has evolved from 
being mainly a theoretical area of Physics to become a 
field supported by observational data of high-precision. 
Recent experiments call upon state of the art technology 
in Astronomy and in Astrophysics to provide detailed in- 
formation about the contents and history of the Universe. 
This has led to the measuring of the parameters that de- 
scribe our Universe with increasing precision. In this 
context, the Nobel Prize of Physics 2011 was attributed 
for the discovery of the accelerating expansion of the Uni- 
verse through observations of distant supernovae ,1,]. In 
fact, after fifteen years of extensive research, we still lack 
a fundamental understanding of this late-time cosmic ac- 
celeration, and its resolution will be one of the major 
objectives in cosmology during the next decade with sev- 
eral surveys and experiments to address the nature of 
the cosmic speed-up. Recently, the released Planck data 
of the 2.7 full sky survey [3| have also shown a number 
of intriguing features, which will require a change in our 
standard view of the Universe. 

The standard model of cosmology has favoured a miss- 
ing stress-energy component, in particular, the dark en- 
ergy models Q, although one may also explore the al- 
ternative viewpoint, namely, through a modified grav- 
ity approach, in particular, using generalizations of the 
Einstein-Hilbert Lagrangian [343 ■ fact, it is funda- 



mental to understand how one may differentiate these 
modified theories of gravity from dark energy models. 
Note that generalizations of the action functional can be 
approached in several ways. For instance, prescriptions 
consist in replacing the linear scalar curvature term in 
the Einstein-Hilbert action by a function of the scalar 
curvature, f{R) gravityM, or by more general scalar in- 
variants of the theory j5| , or still considering a coupling 
between curvature and matter Q, a class of theories of- 
ten termed higher-order gravity theories. Different ap- 
proaches have been analysed in the literature: The metric 
formalism f4j, which consists in varying the action with 
respect to g^''. The Palatini formalism [1, 0|, where the 
metric and the connections are treated as separate vari- 
ables; and the metric-affine formalism, where the matter 
part of the action now depends and is varied with respect 
to the connection [9]. 

A novel approach to modified theories of gravity was 
recently proposed, consisting of adding to the metric 
Einstein-Hilbert Lagrangian an f{TZ) term constructed 
a la Palatini [lo| . It was shown that using the respective 
dynamically equivalent scalar-tensor representation, the 
theory can pass Solar System observational constraints 
even if the scalar field is very light. This implies the ex- 
istence of a long-range scalar field, which is able to mod- 
ify the cosmological [11] and galactic dynamics [l2|, but 
leaves the Solar System unaffected. Static and spher- 
ically symmetric solutions, in particular, wormhole so- 
lutions, were analysed in 13|. A generalization of the 
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hybrid metric-Palatini gravitational theory was also pre- 
sented ^ii], where the gravitational action is taken to 
depend on a general function of both the metric and 
Palatini curvature scalars. The dynamical equivalence 
with a non-minimally coupled bi-scalar field gravitational 
theory was proved and the evolution of cosmological so- 
lutions was studied using dynamical systems techniques. 
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In fTs*! the ghost analysis for these theories was recently 
performed showing that this biscalar generalization al- 
ways presents either ghosts or tachyons. In fact, new 
types of second, fourth and sixth order derivative gravity 
theories were investigated and the hybrid metric-Palatini 
theory was singled out as a viable class of "hybrid" ex- 
tensions of General Relativity. 

The cosmological applications of this hybrid metric- 
Palatini gravitational theory were also further ex- 
plored [ll|. More specifically, cosmological solutions 
coming from its scalar-tensor representation of hybrid 
metric-Palatini gravity were found and criteria to ob- 
tain the cosmic acceleration were discussed. Further- 
more, several classes of dynamical cosmological solutions, 
depending on the functional form of the effective scalar 
field potential, describing both accelerating and deceler- 
ating Universes were explicitly obtained. The cosmologi- 
cal perturbation equations were also derived and applied 
to uncover the nature of the propagating scalar degree of 
freedom and the signatures these models predict in the 
large-scale structure. 

Relative to the galactic dynamics, the generalized 
virial theorem and certain astrophysical applications in 
the scalar-tensor representation of the hybrid metric- 
Palatini gravity were explored in [T^]. More specifically, 
taking into account the relativistic coUisionless Boltz- 
mann equation, it was shown that the supplementary 
geometric terms in the gravitational field equations pro- 
vide an effective contribution to the gravitational poten- 
tial energy. It was also shown that the total virial mass 
is proportional to the effective mass associated with the 
new terms generated by the effective scalar field, and the 
baryonic mass. This fact demonstrates that the geomet- 
ric origin in the generalized virial theorem may account 
for the well-known virial theorem mass discrepancy in 
clusters of galaxies. Thus, hybrid metric-Palatini grav- 
ity provides an effective alternative to the dark matter 
paradigm of present day cosmology and astrophysics. 

In this work, in the cosmological context, motivated by 
the fact that the Einstein static Universe has always been 
of great interest in various gravitational theories, and 
inspired by the possibility that the universe may have 
started out in an asymptotically Einstein static state 
in the inflationary universe context [l^, we analyse the 
Einstein static Universe and its stability in the hybrid 
metric-Palatini gravitational theory. In fact, the Einstein 
static universe has been analysed in a plethora of con- 
texts, such as: in General Relativity with a non-constant 
pressure p/7j; in brane world models [l^; in Einstein- 
Cartan theory [l^; in loop quantum cosmology 20]; in 
f{R) modified theories of gravity [2T| - [23| : in modified 
Gauss-Bonnet gravity [l^l; in IR modified Hofava grav- 
ity [25]; and in massive gravity [26| . 

This present paper is outlined in the following man- 
ner: In Sec. ini we briefly review the formalism of the 
hybrid metric-Palatini gravitational theory, in particular 
the action and field equations, and in Sec. lIIII present the 
respective modified Friedman equations. In Sec. IIVI we 



consider linear homogeneous perturbations in the context 
of the Einstein static Universe in the above-mentioned 
hybrid metric-Palatini theory, and analyze the respective 
stability regions. In Sec. [Vj we summarize and discuss 
our results. 



II. HYBRID METRIC-PALATINI GRAVITY: 
FORMALISM 

The action for the hybrid metric-Palatini gravity is 
1 

2^ 



5* 



d^xy^[R + fin)]+ S„ 



(1) 



where = SttG and we set c — 1. Sm is the matter ac- 
tion, R is the metric Einstein-Hilbert term, IZ = g^^lZ^y 
is the Palatini curvature. IZ^i, is defined in terms of an 
independent torsion-less connection F^^ as 



V 



f L • (2) 



Varying the action ^ with respect to the metric, one 
obtains the following gravitational field equations 



where the matter stress-energy tensor is defined as 
^ _ 2 6{^£„,) 



'9 ^gt"- 



(3) 



(4) 



The independent connection is compatible with the met- 
ric _F(72.)g^i/, which is conformal to 5^^, and where the 
conformal factor is given by F{Tl) = df{TZ)/dTZ. The 
latter considerations imply that 



1 



' 2F2(7^) 

1 



F{n) 



1 1 



2F{Tl) 



<?^,□F(7^) . (5) 



Note that TZ can be obtained from the trace of the field 
equations ([31), which yields 



F{n)n - 2f{n) -r = k^t. 



(6) 



Introducing an auxiliary field, the hybrid metric- 
Palatini action ([1]) can be turned into a scalar-tensor the- 
ory given by the following action (we refer the reader to 
for more details) 



1 



2k2 



(7) 



Varying this action with respect to the metric, the scalar 
(j) and the connection yields the following field equations 



1 




(8) 
(9) 
(10) 
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respectively. Note that the sohition of Eq. (ITOl) imphes 
that the independent connection is the Levi-Civita con- 



nection of a metric h,, 



Thus we are deahng 



with a bi-metric theory and TZfj,i, and are related by 



1 



and consequently 



(11) 



3 



which can be used in the action ([7]) to eliminate the in- 
dependent connection and obtain the following scalar- 
tensor representation 



1 



{l + (t>)R- 



_3_ 

24> 



(13) 



It is important to note that this action differs fundamen- 
tally from the w — —3/2 Brans-Dicke theory due to the 
coupling of the scalar to the curvature. The trace of 
Eq. dl]) yields ~R~ (j)R + 2V ^ k^T, and using Eq. ([9]), 
takes the following useful form 



2V - (j}V^ = K^T + R 



(14) 



Now substituting Eq. dH) and Eq. ([TTl) in Eq. (0), the 
metric field equation can be written as an effective Ein- 
stein field equation, i.e., G^^ 
five stress-energy tensor is given by 



K^T^^, where the effec- 



rr-icff 



9^^ iV + 2a(, 



^-^d^cbd.cb+^g.^id^f] }.(15) 



The scalar field is governed by the second-order evolu- 
tion equation (we refer the reader to fioj for more details) 



1 

2^^ 



<I>[2V - (1 + <I>)V^] 



^T, (16) 



which is an effective Klein-Gordon equation. Equa- 
tion (|16l) shows that the scalar field is dynamical, con- 
trary to the Palatini case where w — —3/2. Thus, the 
theory is not affected by the microscopic instabilities that 
arise in Palatini models with infrared corrections. 



III. EINSTEIN STATIC UNIVERSE AND 
MODIFIED FRIEDMANN EQUATIONS 

A. Effective Friedmann equations 

The cosmological applications of the hybrid metric- 
Palatini gravitational theory have been extensively ex- 
plored [lO, in the scalar-tensor representation. In 



particular, it was shown that the accelerating expansion 
and, in general, any cosmological behaviour, strictly de- 
pends on the effective scalar field potential. However, it 
is rather important to stress that the scalar field has a 
purely geometric origin and describes further degrees of 
freedom of the gravitational field coming from extended 
theories of gravity. 

Now, consider the Friedmann-Robertson- Walker 
(FRW) line element given by 



(12) ds=-dr + a\t) 



1 - Kr"^ 

where K = -1-1,0,1 corresponds to a closed, flat, and 
open universe, respectively. For this case, the scalar cur- 
vature takes the following form 



R = Q[2H 



H 



(18) 



The effective Friedmann equations can be written in 
terms of the effective energy density and pressure, given 
by 

H = (Poff + Pcff ) + , 



(19) 
(20) 



with the following relationships 

(1 + 0) K^pcff = 
(l-f (/)) K^Pcff = 



-Y^<i? + \y{^) - + K^p„, , (21) 



j,2 



1 



V{<j)) + + 2H(j) + K>„ 



(22) 

The conservation equations for the matter component 
and the scalar field are 



3H</)- 



02 

2(j) 



Prn+iH{p„i+p,n) = 0, (23) 

+ l(l>R-l;(l>V'{(f>) = 0. (24) 



B. Einstein static Universe 

For the Einstein static universe we choose a = uq — 
const, so that H = H ^ 0, and the curvature scalar re- 
duces to i? = GK/gq. Note that we do not specify K at 
this point. We also consider that the matter distribution 
obeys the following linear equation of state pm = wpm- 
In classical General Relativity, the presence of a posi- 
tive cosmological constant allows us to solve the analogue 
equations of and (^0)1 under the staticity assump- 
tions with K = 1. With the above considerations, the 
effective Friedmann equations (fT9)) and (|20t reduce to 

^, (25) 



-TT (Poff + Pcff) 

2 a, 



K 



(26) 
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respectively, which imposes the following condition on 
the distribution of matter 



Pcff + 3pcff = . 



(27) 



As a first approach, let us now assume that (f) = (pQ = 
const, then the condition on the effective matter equation 
implies 



and the modified Klein-Gordon equation 
the following expression 



(28) 



reduces to 



6K 



(29) 



Therefore, we can express pm and uq in terms of and 
the potential V. Substituting back into one of the two 
field equations, we can find an implicit equation which 
fixes the value of the field 0o by 



(30) 



Depending on the functional form of the potential, it may 
be possible to solve for (po explicitly. 

One should also note at this point that A' = — f is 
a possible parameter for an Einstein static universe, in 
the context of the hybrid metric-Palatini gravitational 
theory. Provided V'{(po) < 0, one can find the scale factor 
ao using Eq. (f^ . 



IV. STABILITY ANALYSIS OF THE EINSTEIN 
STATIC UNIVERSE 

It is interesting that the Einstein static universe has 
been revived as the asymptotic origin of an emergent 
universe, more specifically as an inflationary cosmology 
without a singularity T^. Despite the fact the positive 
curvature is negligible at late times, its role is crucial 
during the the early universe. In the latter context, it 
allows these cosmologies to inflate and later reheat to a 
hot big-bang epoch. In fact, these cosmological models 



possess attractive features such as the absence of a sin- 
gularity, of an 'initial time', of the horizon problem, and 
the quantum regime can even be avoided. 

In order to study the stability of the Einstein universe, 
we now consider the following homogeneous perturba- 
tions 



a = ao + Sa{t) , 
(l)^(f>o + S(f>{t) , 



P 



(0) 



(31) 
(32) 

(33) 



where we only consider adiabatic perturbations 6pm {t) = 
wSpm{t)- When perturbing the field equations 
and (1^ . and the evolution equations for the matter (^5)) 
and (j24p . we will only consider linear perturbations and 
neglect all terms of higher order. 

Let us begin with equation (|23p . Perturbing to linear 
order provides 



(0) 

Spin + + w) -^Sa = , 
ao 



(34) 



which can be integrated to yield the following relation- 
ship 



(0) 

Spm = —3(1 -I- w) Sa , 

ao 



(35) 



which can be used to substitute the matter perturbation. 

Note that we have three independent equations for 
the three quantities Sa, Scj) and Spm- Thus, we have 
two more equations to solve, these are equations (j20p 
and We are now linearising these with respect 

to the perturbations and substitute (1351) and the back- 
ground solutions p5|) -(|30 |) into the field equations (|20|) 
and (|24ll . This results in two second order ordinary dif- 
ferential equations. It is useful to write the result in the 
following form 



Sa 



Mil Mi2 
M21 M22 



Sa 



(36) 



where the components of the matrix M are given by 



Ml 
Ml 

M2 

M- 



22 



_ 3y((/)o)(l + w){-l + iw) + 2r (0o)(l + 3u;) 
6(1 -fSw) 

_ /^ -3y(0o)F"(0o)(l + w) + V'{(t)o)[V'{^o) + V"{cpo)]{l + iw) 

V6F'(,^o)3/2(l + 3w) 
^ V{MV'{M^'''{'^ + w)(-l + 3w)[r(0o)(l + M - 3l^(0o)(l + w)] 

^ [r(0o)(l + 'iw) - W{4>o){l + w)][^W{My"{4>o){y + + V'{4>oY{l + 3w)] 

3y'(0o)2(l + 3w)2 



(37) 
(38) 
(39) 
(40) 
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1 






S(f> 










^) 




Sb 


Mil 


M12 







Sa 




VM21 


M22 





0/ 





It is more convenient to work with a first order than 
a second order system of equations. Thus, introducing 
two new variables, Sb = da and dtp = 5(j), Eq. ([36|) can be 
rewritten as 



(41) 



This last equation is of the form v — A.V , and the eigen- 
values of the 4x4 matrix A determine the stability or 
instability properties of the perturbations. Let us denote 
the two eigenvalues of M as Ai and A2, then the four 
eigenvalues of A are simply given by ±y/Xi and ±-\/A2. 
We find 

Al - i (Mh + M22 + V(Mll - Af22)' + 4M12M21) , 

(42) 



A2 = - (Afii + M22 - ViMii - M22f + 4M12M21 

(43) 

In order for stable perturbation to exist, we need to re- 
quire that both eigenvalues are purely real, and addition- 
ally Ai < and A2 < 0. The reality condition entails 



(Afii - M22) + 4Afi2A/2i > , 
while the negativity condition implies 



Mil +M22 < - V(Mii - M22Y + 4Ari2M2i , 
(Mil + M22)^ < (Mil - M22)^ -I- 4Ari2Af2i . 

The latter conditions can be rewritten as 

{Mil + M22f - (Mil - M22f < 4Afi2Af21 , 

4AfiiM22 < 4M12M21 , 
M11M22 - M12M21 < , 
det M < . 



(44) 



(45) 
(46) 



(47) 
(48) 
(49) 
(50) 



Therefore we have the necessary condition detM < 0. 

Taking both conditions together, namely the reality 
and the negativity conditions, we require 

(A/11 - Af22)^ + 4Afi2M2i > (51) 
A/11 + A/22 < - V(Mii - A/22)2 + 4A/12A/21 . (52) 



As our expression for the components of M are 
quite involved, it is expected that those inequali- 
ties are very complicated when expressed in terms of 
V{(l)o), V'{(j)o), V"{(j)o) and w. However, we can indicate 
graphically the parameter space where the conditions are 
satisfied. Firstly, we note that we can always rescale 
the values of the potential at (/fp, therefore we will set 
V{(j)o) = 1. For a fixed value of w, the two conditions 
now depend on the two values V'{4>o) and V"{(f>o)- We 
plot the relevant regions in the Figures [T]-[31 



In Figures HHSl the shaded regions represent the param- 
eter space where homogeneous perturbations of the Ein- 
stein static universe are stable. In each figure the equa- 
tion of state parameter w is fixed, we chose the six rep- 
resentative values w — {1/2, 1/5, 0, —1/5, —3/5, —4/5}, 
respectively. Note that a wide variety of stability regions 
exist for a range of positive and negative values of V'{(I)q). 
One could have attempted a three dimensional plot with 
axes V {(f>o),V" {4'o),w, however, this gives a very com- 
plicated figure of little value. We should emphasise that 
these stability figures have to be seen in context with the 
background solutions (P5)) - ([5(7|) . For instance, we require 
a positive effective energy density in (P5|) and a positive 
value for ao in ([29)) . both of which are independent of 
the stability of the perturbations. We should also note 
that once a specific potential has been chosen, the terms 
F'((/)o) and F"((/)o) are not necessarily independent any- 
more and care has to be taken when interpreting the 
figures. 



V. SUMMARY AND DISCUSSION 

In this work, we analysed the stability of the Ein- 
stein static Universe by considering linear homogeneous 
perturbations in the respective dynamically equivalent 
scalar-tensor representation of hybrid metric-Palatini 
gravity. Considering a linear equation of state parameter 
for the matter distribution, the stability regions of the 
Einstein static universe were parametrized by the first 
and second derivatives of the scalar potential, and it was 
shown that a large class of stable solutions exists in the 
respective parameter space. The stability regions were 
considered by requiring positive values for the effective 
energy density and for the scale factor ag, both of which 
are independent of the stability of the perturbations. 

It is interesting to discuss a few specific choices for the 
potential V{(l)) and search for stable static solutions. For 
instance, consider the exponential potential given by 



F(0) = Foe-^^ 



(53) 



where Vq and A are two positive constants. In an open 
universe (K = — 1) a static solution is given by 



ao 



Pq 




exp 



A 3(A l)u; -I- 3 



exp 



A + 



6w + 2 
2 



A + 3(A + l)u; + 3 
A(l + 3w) 



3m; -h 1 



(54) 
(55) 
(56) 



where w > —1/3 in order to have po > 0. This is quite 
interesting as the effects of hybrid metric-Palatini gravity 
are sufficiently strong to allow for a static solution in the 
open case. Unfortunately this solution is always unstable 
for w > —1/3. This is in line with the point made above. 



V'(il>o) 



FIG. 1: Left: Stability regions for w = 1/2. Right: Stability regions for w — 1/5. The shaded regions of the figures represent 
the parameter space where homogeneous perturbations of the Einstein static universe are stable. In both cases one cannot find 
stable solutions unless V^'(<j!>o) < 0. 
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FIG. 2: Left: Stability regions for w = 0. Right: Stability regions for w — —1/5. As in the previous figure, the shaded regions 
of the figures represent the parameter space where homogeneous perturbations of the Einstein static universe are stable, and 
in both cases the stability regions exist for V'{4>o) < 0. 



namely that the stability regions have to be interpreted 
using specific background solutions. 

Let us now set if = 1. Recall that in General Rela- 
tivity the Einstein static universe only exists ii K = 1. 
Kowever, for this case, there are no real and positive solu- 
tions for oo in this closed universe with the exponential 
potential ([53l) . This follows from (|29p . the exponential 
potential has negative derivative and for positive K the 
scale factor oq cannot be real. This indicates that one 



should consider a different potential. 

As a second specific potential, consider the quadratic 
potential given by 

V{ct>) = \m'<i? , (57) 
where to is a positive constant scalar field mass. A static 
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FIG. 3: Left: Stability regions for w = —3/5. Right: Stability regions for w — —4/5. The shaded regions of the figures 
represent the parameter space where homogeneous perturbations of the Einstein static universe are stable. The equation of 
state parameter now satisfies w < —1/3 and one can find stable solutions for a wide variety of positive values of V'{(j}o)- 



solution for this potential is given by 



ao 



Po = 



3k{l - 3w) 
Y A(6w + 2) 

4:A{3w + l) 

(1 — Sw)^ ' 

6w + 2 
1 - 3w ' 



(58) 

(59) 
(60) 



As in the previous case, the background solution will 
place constraints on the permissible parameter range. In 
an open universe {K = —1) we must require w > 1/3 
in order to have ap and po both real and positive, while 
in a closed universe {K — 1) we need —1/3 < w < 1/3. 



The potential satisfies V"{4>) 



and is thus strictly 



positive. Combining this with the condition on w, we 
conclude that the open universe solution is always unsta- 



ble. On the other hand, the closed universe is stable for 
every allowed value of w. 

Despite concentrating the analysis in this work to 
the stability of the Einstein static Universe by consider- 
ing linear homogeneous perturbations, which has yielded 
interesting results and has presented a richer stabil- 
ity/instability structure that in General Relativity, it is 
of interest to extend our results to inhomogeneous per- 
turbations. Work along these lines is presently underway. 
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